arXiv: 1506.05501 v3 [astro-ph.EP] 25 Nov 2015 


Draft version November 26, 2015 

Preprint typeset using style emulateapj v. 5/2/11 


ATMOSPHERIC CHEMISTRY FOR ASTROPHYSICISTS: 

A SELF-CONSISTENT FORMALISM AND ANALYTICAL SOLUTIONS FOR ARBITRARY C/O 

Kevin Heng', James R. Lyons^, and Shang-Min Tsai* 

Draft version November 26, 2015 

ABSTRACT 

We present a self-consistent formalism for computing and understanding the atmospheric chemistry of exo¬ 
planets from the viewpoint of an astrophysicist. Starting from the hrst law of thermodynamics, we demonstrate 
that the van’t Hoff equation (which describes the equilibrium constant), Arrhenius equation (which describes 
the rate coefficients) and procedures associated with the Gibbs free energy (minimisation, rescaling) have a 
common physical and mathematical origin. We address an ambiguity associated with the equilibrium constant, 
which is used to relate the forward and reverse rate coefficients, and restate its two definitions. By necessity, 
one of the equilibrium constants must be dimensionless and equate to an exponential function involving the 
Gibbs free energy, while the other is a ratio of rate coefficients and must therefore possess physical units. We 
demonstrate that the Arrhenius equation takes on a functional form that is more general than previously stated 
without recourse to tagging on ad hoc functional forms. Finally, we derive analytical models of chemical 
systems, in equilibrium, with carbon, hydrogen and oxygen. We include acetylene and are able to reproduce 
several key trends, versus temperature and carbon-to-oxygen ratio, published in the literature. The rich variety 
of behavior that mixing ratios exhibit as a function of the carbon-to-oxygen ratio is merely the outcome of 
stoichiometric book-keeping and not the direct consequence of temperature or pressure variations. 

Subject headings: planets and satellites: atmospheres - methods: analytical 


1. INTRODUCTION 
1 . 1 . Preamble 

Understanding chemistry is indispensible to deciphering 
the abundances of atomic and molecular species present 
in an exoplanetary atmosphere. Despite its somewhat late 
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atmospheric chemistry has a long and rich history in 
the Earth and planetary science s and t h e study of brown 
dwarfs (e.g., [Prinn & Barshavl |1977t iBarshav & Lewi: 
WtI f^dlCT&Yungl II98Tr lFeglev & Fodders! I~ 


Fodders & FeglevT 2()02l : TCiesla & CharnlevI I 2 OO 6 I) . Yet, 


a first-principles, self-consistent formalism that unihes all 
of the quantities and terminology in a form that is useful 
for astrophysicists is missing from the literature. For exam¬ 
ple, there is more than one definition of the “equilibrium 
constant”. 

Within the same framework, we demonstrate that the van’t 
Hoff equation (which describes the dimensionless form of 
the equilibrium constant), the Arrhenius equation (which de¬ 
scribes the rate coefficients) and procedures associated with 
the Gibbs free energy (minimisation and scaling) all origi¬ 
nate from the hrst law of thermodynamics. The foundations 
of atmospheric chemistry are built upon statistical mechan¬ 
ics, since the hrst law derives from it. To demonstrate the 
usefulness of our formalism, we use it to compute analyt¬ 
ical solutions of chemical systems with pure hydrogen and 
with carbon, oxygen and hydrogen (gas ph ase only); we show 
that these solutions generalise the work of iBurrows & Shaip 
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(I1999^ and correctly reproduce all of the expected trends. 


1 . 2 . Survey of monographs 

The novelty of the present study is not in the individual for¬ 
mulae stated, which are mostly previously known; with one 
exception, we certainly do not claim to be the hrst to derive 
these formulae. Rather, it is in the way these results are de¬ 
rived and unihed under a common, self-consistent, mathemat¬ 
ical formalism that is accessible and palatable to astrophysi¬ 
cists (rather than to chemists). We will now demonstrate this 
claim of novelty by surveying several textbooks in chemistry. 

Specihcally, equations (|5]l, dTji, (fTOl i. (fOl) . (fTTl) and (l 22 l l 
are commonly stated in textbooks. Our intention is to 
weave a common mathematical th read be t ween t hem. We 
have s ur veyed the monographs o f ISlate d (ll^39|l. riohnstonl 


(I1966^■ I van Zegg eren & StorevI (Il97(]h. Moord (1^72), 


lEise^erg & CrothersI (Il979h . ISmith &Missenl (1982 ) 


fSteinfeld . Francisco & Hasel (1198^ . lAtkins & de Paula 
l2006h. iKlotz & Rosenberl (l 20 m IPeVod (1201 5h and 
Glassman, Yetter & Glumad (120151) and verihed that, while 
each lists some subset of these formulae, none of them derive 
and unify all of these formulae in the manner of the present 
study. None of these monographs derive the generalised form 
of the Arrhenius equatio n that we present in equat ion (l22l l or 
the generalisation of th e [Burrows & Sharpl (119991) analytical 
solutions we present in 114.21 


2. GENERAL SETUP 

In general terms, we consider a chemical reaction involv¬ 
ing a pair of reactants (Xi and X 2 ), which produces a pair of 
products (Zi and Z 2 ), 

oiXi -f 02X2 ^ biZji -f 62Z2, (1) 

where ai, 02 , 61 and 62 are the stoichiometric coefficients. 
The reactants and products may be atoms or molecules of 
arbitrary stoichiometry. The forward and reverse reactions 
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are described by the rate coefficients fcf and fcr, respectively. 
Whenever we discuss something in general terms, it will al¬ 
ways be with reference to the preceding chemical reaction. 


3. EQUILIBRIUM CHEMISTRY: GIBBS FREE ENERGY, 
EQUILIBRIUM CONSTANT(S), VAN’T HOFF’S EQUATION AND 
ARRHENIUS’S EQUATION 

3.1. The Gibbs free energy 

If we denote the specific internal energy by U, the tem¬ 
perature by T and the specific entropy by S, then a reason¬ 
able guess for the excess energy associ ated with a chemical 
reaction is the Helmholtz free en e rgy (ISlateHll939l: Moorel 

ins & de Paula 120061 


wn lEisenberg & Crothersll l979I lAtkii 
S wendsenll20 1 2i: iDe Voell20 1 sIl 


F = U-TS. 


( 2 ) 


It turns out that this quantity is not general enough be¬ 
cause it does not consider the work done on the sys¬ 
tem. The general quantity is known as the Gibbs free 
energy (ISlated II939I: Ivan Zeggeren & Store\^ Il970t iMoorel 


1972t lEisenberg & Crothers 


Swendsenll20I2[lDeVoell20T¥ 


19791 lAtkins & de PauIail2(M 


G = F + PV, 


(3) 


where P is the pressure, V = l/p is the specific volume and 
p is the mass density. 

The Gibbs free energy plays a role analogous to the La- 
grangian of classical mechanics, which is the difference be¬ 
tween the kinetic and potential energies of a system. Instead 
of solving Newton’s equation directly, one may minimise the 
Lagrangian, a technique known as the principle of least ac¬ 
tion. Gibbs free energy minimisation and chemical kinetics 
are the chemical analogues to these two techniques. 

In a chemically-active system, the number of particles of 
each species is generally not a conserved quantity. If we de¬ 
note the number of particles associated with the j-th species 
of the system by Nj , then the first law of thermodynamics 

; JacobsonI 


needs to be modified ( van Zesaeren &_Storev 2970t Jacobs 
l2005l:f^endsenll2012[lGlassman. Yetter & Glumad 201 fill . 


TdS = dU + P dV - '^CjdNj , (4) 

3 

where Cj is the chemical potential associated with each 
species. The sum is performed over all of the species in the 
system. 

By usin g the definition of G and the product rule, one may 
show t hat (ISmith & Missenlll982t iKlotz & Rosenberg 120081: 

iDeVoellinsli 


dG = VdP - SdT + Y, Cj dNj. (5) 

3 


We will now proceed to show that a variety of useful quanti¬ 
ties originate from this equation, which is essentially still the 
first law of thermodynamics. 


3.2. Gibbs free energy minimisation 

Generally, the entropy of a system increases according 
to the second law of thermodynamics; at constant tempera¬ 
ture and pressure, its Gibbs free energy generally decreases 
and seeks a minimum. If one is interested in solving for 
the chemical equilibrium of a network of reactions, then 
one needs to minimise the Gibbs free energy of the system 


(Ivan Zeggeren & Storevl[T970h . Within the context of our for¬ 
malism, we will now elucidate the exact expressions involved 
in this minimisation. Generally, we have Cj = Cj (T, P) and 
equation (|5]l cannot be straightforwardly integrated. However, 
at a constant temperature and pressure—which is the typical 
circumstance under which one performs Gibbs free energy 
minimisation—equation (|5]l reduces to 

dG ^Y^3^^P ( 6 ) 

3 

The integration can be p e rformed triviall’v to yield 
(Ivan Zeggeren & Storeyl 119701: lEisenberg & CrothpsI 119791 : 

ISmith & Missenl|1982l:lGlassman. Yetter & Glumacll201^ 

G = (7) 

3 

It is not uncommon to see G being defined as the product of 
the chemical potential and the number of particles of a given 
species. Strictly speaking, it is not a definition—rather, it is 
the expression for G in the isothermal and isobaric limit. 

Equation (|7]i is the quantity we need to minimise, but we 
need additional equations to close the system. In the absence 
of nuclear reactions, this arises naturally from the notion that 
the elemental building blocks of molecules cannot be created 
or destroyed. Thus, the number of carbon, hydrogen, oxygen, 
etc, atoms in a system is invariant between the reactants and 
the products, whether they exist in their atomic form or are 
sequestered in molecules. Mathe matically, this set of b ook- 
keeping equations takes the form (ISmith & Misseniri982f) . 

YA^,N,=N'. (8) 

3 

The matrix Aij states the number of atoms of species i present 
in the molecular species j. The number of atoms of species i 
is denoted by Nf 


3.3. The equilibrium constant: more than one definition 


A persistent source of confusion exists in the literature re¬ 
garding the definition of the equilibrium constant. Seyeral 
references list it as being composed of a series of partial pres¬ 
sures associated with the reactants and products and equates it 
to an exponential term i n yolying the Gibbs free energy (e.g.. 
Burrows & Sharnl Il999t iJacobsonI 120051: IVisscher & Moses 
201TI iKopparapu, Kasting & Zahnlel 1201 3 iLine & Yung 

ISj, 


K, 


eq,literature — 


pai pa2 • 
Xi X2 


(9) 


Taken at face yalue, the partial pressure has physical units. 
Since the stoichiometric coefficients of the reactants ( 01 - 1 - 02 ) 
and the products (61 + 62 ) are generally unequal, TTeq,literature 
must generally haye physical units and cannot be equated to 
an exponential term (which is by definition dimensionless), 
unless the partial pressures haye somehow been normalised. 
It is not always explicitly explained that this normalisation 
has been performed. Seyeral monogra phs haye preyiously 
mentioned this normalisation procedure (|Fermi|| 1936t IMoorel 
119721 ISmith & Missed IT982I : lAtkins & de Paulal I2006h . but 
we will now proyide a deriyation that is consistent with the 
rest of our formalism. 

To deriye the equilibrium constant, we return to equation 
(|5]l and consider it in the limit of dT = 0 and dNj = 0. If 
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we invoke the ideal gas law (P = pTZT), t hen we obtain (e.g., 
lEisenberg & Crotherslll979l:lDeVoell2015h 


G = Go +7^Tln 


( 10 ) 


where Pq is a reference pressure, TZ is the specific gas con¬ 
stant and Go = G{Po,T). The preceding equation is useful 
for scaling the Gibbs free energy to other pressures given its 
value at a reference pressure—it is exactly the equation one 
has to use when extracting G for different values of P from 
thermodynamic databases, which typically tabulate values of 
Go. 

We may use the preceding expression to combine 
the Gibbs free energy of the reactants and pr oducts, 
weigh t ed by their stoichiometr i c coefficients ( Slater! 

lEisenberg & Crother'^ l^th & Missed iT^ : 

Klotz & Rosenbergl2008l:lGlassman. Yetter & Glumacl2015h 


AGo - AG = -7^^1niT, 


eqi 


where we have defined 


AG = 6iGzi + & 2 GZ 2 — oiGxi — a2Gx2, 

AGo = &iGzi,o + (' 2 Gz 2 ,o — aiGxi,o — a 2 Gx 2 ,o. 


( 11 ) 


( 12 ) 


What is interesting is that this first-principles approach nat- 
urally yields the definitio n for the equilibrium constant (e.g., 
iKlotz & Rosenber3l2008l) . 


f)2 


_ jp^jPoriPzjPo) 

- {PxjPor {PxjPor 


(13) 


Notice that this equilibrium constant is naturally dim ension¬ 
less; i ts deriyation is similar t o the ones giyen in iMoorel 
iSmith & Missed (fl9^ . lAtkins & de Paulal (120061) 
and lDeVoel (120151) . who obtained it in terms of chemical po¬ 
tentials and actiyities. The factors of Pq appear without being 
inserted in an ad hoc manner. 

Physically, the system adjusts itself until it reaches 
chemical equilibrium, which occurs when AG = 0. Let 
the reference state, characterised by Pq, not be in equi¬ 
librium, such that AGo ^ 0. If one is referring to a 
molecule, then AGo is the energy needed to construct it 
from its constituent atoms—it is the Gibbs free energy of 
formation. If one is referring to mixtures of molecules, 
then AGo is the difference in the Gibbs free energies of 
formation between the reactants and products. Equation 
TTb na t urally y ie lds the re la tionsh ip between AGo and Kf.„ 


19^ 
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119821 
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1989, 

Atkins & de Paulal 120061; iKlotz & RosenbergI 120081 iDeVoe 


Ken = exp - 


AGo 

7^T 


(14) 


Eor example, while equations (A l) and (A6) of 
[Burrows & Sharp! (!1999!). equation (7) of IVissc her & Mose^ 
(1201 ll) , equation (5) ofjKopparapu , Kas t ing & Zahnle! (12012!) 
and equation (2) of iLine & Yung! (1^1^ do not explicitly 
mention the factors of Pq needed to render ATeq dimension¬ 
less, it is common practice to omit these reference-pressure 
terms in standard treatments of chemical equilibria. 

We now seek another possible definition of the equilibrium 
constant. Let the number density be generally represented 


by n; self-explanatory subscripts relate it to the appropriate 
reactant or product. In chemical equilibrium, we expect the 
forward and r eyerse rate coefficients to be related as follows 
(IMoorel! 1972!) . 

Wi»^X2 = (15) 

A plausible, alternatiye definition for the equilibrium constant 
is i!Johnstonl!l966l : !Moore![l972t[^infeld, Erancisco & Hasel 

!T^ 

( 16 ) 

Note that since kf and kr generally do not possess the same 
physical units, AT'q is expected to be dimensional. 

We may relate our two defin itions of the equilibrium con¬ 
stant (lAtkins & de Paul^!2006!) . 

AT'q = Keq p&i+«'2-ai-a2 ^ (27) 

with kp, being the Boltzmann con stant. As has been pointed 
out by IVisscher & Moses! (1201 Ih . the “pressure correction 
term” (which is really a temperature correction term) is some¬ 
times missed by other workers. It is less well- known that thi s 
correction term has already been elucidated by lEermi! ([193^. 
It yanishes when ai + 02 = bi + 62 and we haye ATeq = K'eq 
(which is commonly, but not always, true). 

In other words, K' is used to “reyerse” the forward rate 

“4 

coefficients, but it is iTeq that relates it to the Gibbs free en¬ 
ergy. Equation (ITtI i relates them properly. In the literature, 
what we haye defined as iTeq and AT'q are often denoted, re- 
spect iyely, by Kp and Kgq instead (e.g.. IVisscher & Moses! 
120111). although such an approach is not uniyersally adopted 
(e.g.. !Line & Yungl!2013!) . 

3.4. The van’t Hoff equation 

If we differentiate equation (ITTI l with re spect to the tem¬ 
perature, we obtain th e yan’t Hoff equation (ISmith & Missed 

EHmiSEioBIlQQl), 

d (In Keg) _ AGo 

dT ^ 

Note that we are allowed to go from equation (ITTI l to ( fTSl ) only 
because we haye constructed AGq to be isothermal. 

In most incarnations of the yan’t Hoff equation, it is the 
change i n enthalpy, rather than the Gibbs free energy, which 
is stated (lMoorell972HSmith & Missenll982HJacobsonl2005l : 

lDeVoell2015l) . If the system is isoth ermic and ad iabatic, then 

these two statements are equiyalent (lSlateilll93^ . 

3.5. The Arrhenius equation: rate coefficients and activation 

energies 

We next deriye the expressions for the rate coefficients. At 
this point, we need to inyoke the notion of the actiyation en¬ 
ergy, which is the energy barrier associated with a forward or 
reyerse reaction. One may think of the reactants and products 
as being two different stable states residing at different en¬ 
ergy leyels. To transition from one state to the other requires 
that one surmounts an energy barrier, which is the actiyation 
energy (Eigure[T]i. The barrier of the actiyation energy orig¬ 
inates from the need to oyercome bond strengths and the re¬ 
quirement that the reactants haye specific orientations during 
a collision. Eor a single reaction, the difference between the 




































































































4 



Fig. 1.— Schematic depicting the relationship between the activation ener¬ 
gies and the change in enthalpy. 


activation energies of the forward and reverse reactions is the 
change in the enthalpy, 

AGo = Ef-E,-TASo, (19) 


where E{ and E,- are the activation energies associated with 
the forward and reverse reactions, respectively, and ASq is 
the change in entropy at the reference pressure. The preceding 
expression allows us to cast the adjectives “exothermic” and 
“endothermic” in more precise, mathematical terms. If the 
activation energy of the forward reaction exceeds that of the 
reverse one, then one needs to inject energy into the system 
for it to proceed, i.e., Ef — E^ > 0. One refers to this as 
an endothermic reaction. Reactions with E{ — E^ < 0 are 
exothermic. 

By combining the expressions for and iT'q, we obtain 


In fcf — In fcr = — 


E{ — E^ 
TIT 


ASo 

n 


+ (ai + 02 — 6i — 62) In 


k^T 

Po 


( 20 ) 


The symmetries inherent in the preceding equation suggests 
that it may have been constructed from two independent gov - 
erning equations for the rate coefficients (IUpadhvavll2006h . 
Mathematically, “splitting” this equation is a degenerate en¬ 
deavor and is not rigorous. To persist in this endeavor, we 
have to appeal to physics. First, we expect that the governing 
equations for k{ and fcr must enjoy a large degree of symmetry 
between them. Second, we expect k{ and k^ to be associated 
with El and E,;, respectively. Thus, a plausible guess is that 
the preceding equation originated from the difference between 
these two equations. 


In fcf 
In fcr 


El 

-^TcrlnT- 


TIT 


+ Cr In T + 


c'lASo 

7^ 

c'A^o 

n 


Cf 

Cr 


( 21 ) 


The coefficients Cf and Cr cannot be stated uniquely. For ex¬ 
ample, we can have ci = ai + 02 and Cr = -f ^2; we may 
also have Cf = —bi — 62 and Cr = —ai — 02. This mathe¬ 
matical freedom implies that Cf and Cr may take on a range of 
values and may be positive or negative. 


Finally, we end up with the Arrhenius equations, 

k, = A,T« 

where we necessarily have 

Cf — Cr = Oi -I- 02 — 61 — 62 , 

4 - 4 = 1 , 

c" - c" = (oi -f 02 - 61 - 62) In (. 


( 22 ) 


(23) 


The pre-exponential factor Ai absorbs terms associated with 
Cf, c" and ASo', its counterpart, Ar, does the same for c', 
c" and ASq. Absorbing the ent ropy into the pre-expone ntial 
factors was previously noted by I Yung & DeMoi^ (Il999h . but 
our derivation is more general as it involves cj and c'. We note 
that one may also use the van’t Hoff equation as a starting 
point for the derivation. 

Traditionally, derivations or statements of the Arrhe¬ 
nius equation include only the exponentia l term i nvolv¬ 
ing the activation energy ( Johnstoril Il966t iMoore 
IStdnfeld.fiancisco & Has^ 198^ I Yung & DeMore 


1972; 


1999; 


iJacobsonfl^OOSnAtkins & de Pauld 1200^ . They omit the 
powe r-law terms and tag them on, after the fact te.g.. l.Tacobsonl 
l2005h . partially as a means of using them as fitting func¬ 
tions for experimental data. Our deriyation demonstrates that 
there is a sound basis to including these terms. Thus, the Ar¬ 
rhenius equations attain a status that is eleyated aboye that 
of mere ad hoc htting functions. Typically, Cf = 0 and 
Cr = 0 suffice for low temperatures; “non-Arrhenius” behay- 
ior, w here ci A 0 and Cr A 0, is import ant at high tempera¬ 
tures dOlassman, Yetter & Glumacll2015h . 

Kinetic theory states that the rate coefficient is giyen by 
(cTcoiffirei). where CTcoU is the cross section for collisions be¬ 
tween the reactants and Urei is the relatiye yelocity between 
them. If CTcoii is independent of the relatiye yelocity, then 


(ccoiffirei) c)c dOlassman, Yetter & Glumacll2015h . If 
(Tcoii depends on the relatiye yelocity, then more general 
power-law dependences on T are possible. 

The Arrhenius equations do not account for three-body 
reactions. When the number density of the third body is 
low, the reaction rate is linearly proportional to it. As it in¬ 
creases, a point is reached where the reaction rate saturates 
to a limiting yalue. Fitting functions for implementing this 
saturation effect haye pre yiously been giyen by, for example, 
IVisscher & MosesI (1201 Ih . 


4. ANALYTICAL MODELS OF ATMOSPHERIC CHEMISTRY 

The formalism and concepts we haye established may be 
highlighted yia a set of analytical models. 


4.1. Pure hydrogen 

For completeness and as the simplest example, we con¬ 
sider a system consisting purely of hydrogen in its atomic and 
molecular forms. 


2H + M H2 + M, (24) 

where M is a third body of arbitrary stoichiometry. 
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K' = X'q^totai, we may solve for the (normalised) number 
density of atomic hydrogen (iGail & Sedlmavrll2014l) . 


, ^ nn + + 

= - = - J-fy. - 

^total 4rV 


(28) 


The precedi ng expression is similar, but not identical, to that 
presented in lBarshav & Lewlsl (Il978h . 

Since K' oc K'^^P/T, one may argue that increasing K' 
values correspond to decreasing temperatures. Figure|2]shows 
the curves of tt-h and fiH 2 = ^H 2 /^t-totai- As expected, molec¬ 
ular hydrogen prevails at low temperatures. So far, our toy 
model does not allow us to define what “low” is, as we have 
not related K' to T and P. AppendixlBllists the Gibbs free en¬ 
ergies used to make this conversion. In Figure |2] we include 
a separate set of calculations where nn and are shown 
as functions of temperature and pressure. At T < 3000 K, 
hydrogen exists predominantly in its molecular form. 



Fig. 2.— Abundances of atomic and molecular hydrogen, normalised by 
the total number density, as a function of the normalised equilibrium constant 
(top panel) and temperature (bottom panel). In the top panel, K' is a proxy 
for the temperature; larger K' values correspond to lower temperatures. 


Using our formalism for the evolution equations (see Ap- 
pendix| 3 i, we may write down 


1 dn-H. 

2 dt 
i9nH2 

dt 


= + nH 2 ^Mfcr, 


= n^nyiki — riYi^UMkr- 


(25) 


If we add these equations and perform the integration, we end 
up with 

riH -f 2 nH 2 = fftotai- (26) 


This is already a demonstration that the correction factor 
(the reciprocal of the stoichiometric coefficient) i s essential 
(jjohnstonlighhtr^einfeld, Francisco & Hasei:1989l) . if one de¬ 
sires to get the book-keeping between the hydrogen atoms and 
molecules correct. Here, ntotai is the total number of particles 
in the system. 

In chemical equilibrium, the (dimensional) equilibrium 
constant of the reaction is 


T^' _ '^H2 


(27) 


If we plug this expression back into equation (l26T l and define 


4.2. Carbon, hydrogen and oxygen: methane, water, carbon 
monoxide and acetylene 

Inspired by the work of lBurrows & Sha^ ([199^, we seek 
to generalise our toy model of a system with pure hydrogen 
to one that contains carbon (C), hydrogen (H) and oxygen 
(O), albeit only in gaseous form, and any carbon-to-oxygen 
ratio (C/O). We wish to compute the relative abundances of 
the resulting molecules: methane (CH 4 ), water (H 2 O), car¬ 
bon monoxide (CO) and acetylene (C 2 H 2 ). 

We consider the reaction of methane w ith water to form 
carbon monoxide and molecul ar hydroge n dBurro ws & Sht^ 

rr^lTodders & FeglevlUOPTllM OSes et akboilh, 

CH 4 + H 2 O ^ CO -f 3 H 2 . (29) 

The fo rmulae presented in the appendix of lBurrows & Sha^ 
(119991) consider only this reaction and thus are unable to rep¬ 
resent carbon-rich atmospheres, where a variety of hydro¬ 
carbons are present at high temperatures (iLodders & FeglevI 
l2002tlMadhusudhanll2012i:IVenot et al.ll2015h . If these hydro- 

carbons are excluded, then one gets the spurious result that 
methane is always the dominant carbon carrier at high temper¬ 
atures and in carbon-rich situations. Our desire for an analyti¬ 
cal model does not allow us to include all of the hydrocarbons 
that are expected to form. Instead, we assume that acetylene 
is the do minant hydrocarbon and include it yia the f ollowing 
reaction (ILodders & Fegle\^l2002tlMoses et ani201 ll) . 

2 CH 4 t^C 2 H 2 -f 3 H 2 . (30) 

If nitrogen is pr esent, we expect hyd rogen cyanide (HCN) 
to form as well (iMadhusudhanluOl^ . but in the interest of 
algebraic tractabil ity we will not include it. Furthermore, 
I Venot et ^ (l20 1 5l) haye shown using calculations of chemical 
kinetics that acetylene and hydrogen cyanide are the dominant 
hydrocarbons in carbon-rich atmospheres. 

In reality, both reactions are net reactions that consist of 
large networks of indiyidual reactions, some of which pro¬ 
duce transient species en route to the products. We assume 
that hydrogen exists mostly in its molecular form, such that 
the partial pressure of H 2 is, to a good approximation, the to¬ 
tal pressure (P) of the system. Atomic hydrogen is expected 
to introduce only a small correction to P = n-n^^k-BT. This 
simplification essentially remoyes the need for an additional 
equilibrium constant to account for the atomic to molecular 
transiti on (a nd yice yersa) of hydrogen, as was described in 
Section HTT] 
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Fig. 3.— Mixing ratios of methane, water, carbon monoxide and acetylene. 
Again, K' is a proxy for the temperature, but larger K' values correspond 
to higher temperatures. We have set no = 5 x 10““^ as th is is the approx - 
imate value of the Sun’s photospheric oxygen abundance jLodderd 1 2Q03h . 
Top panel: mixing ratios as a function of K'. The “Solai*” and “carbon-rich” 
cases correspond to fiQlhQ = 0.5 and uqIuq = 2, respectively. Bottom 
panel: mixing ratios as a function of the carbon-to-oxygen ratio (nc/f^o)- 
The “cold” and “hot” cases correspond to K' = 10 and K' = 10^, respec¬ 
tively. For illustration, we have set K'^/K' = 10“^. 


Fig. 4.— Same as Figure[3] but with the equilibrium constants being related 
to temperature and pressure via the Gibbs free energies taken from the JANAF 
database. For illustration, we set P = 1 bar. Here, the “cold” and “hot” cases 
correspond to T = 950 K and T = 1500 K, respectively. For the bottom 
panel, note that acetylene has a mixing ratio below 10“^^ (the lower limit of 
the vertical axis) for the cold case. 


The dimensional equilibrium constant of the reaction de¬ 
scribed in equation ( |29] t is 


K^ = 


ncon'^ 




^CHi^HaO JT-CHi^HaO 

while that of the reaction in equation (l30l l is 


K' — 

^eq,2 — 


’^CaHafT-Ha _ ^CaHa^Ha 


''CH 4 


''CH 4 


(31) 


(32) 


Analogous to the case study of pure hydrogen, we have de¬ 
fined 


K' 

K'= 


K' 

_ T^f ^ ^eq,2 

rl2 tl2 


(33) 


as a proxy for the temperature. An important limitation of 
our model is the difficulty with relating K' and K'^, because 
this requires us to explicitly state the functional forms of the 
change in Gibbs free energies of the two reactions. We first 
make the simplest assumption: that K 2 IK' is a constant; we 
will discuss the implications of this assumption later. While 
we could certainly specif y the temperature d e pende nce of K' 
and K 2 —which is what iBurrows & Sharnl (119991) did—we 
initially choose not to so as to understand what such a sim¬ 
ple model would teach us. We will see shortly that the sim¬ 
plicity yields an important insight, which is that we recover 
most of the qualitative trends simply by using the equilibrium 
constants as proxies for the temperature. 

The conservation of particles, as described in equation ([8]), 
states that 


but we note that since K' oc iT'qT^/P^, we expect K' to 
increase with temperature, opposite from the trend associ¬ 
ated with the pure-hydrogen system. We will again use K' 


n-CHi + nco + SncaHs = nc, 

JT-HaO + nco = no, 

4ncH4 + SriHaO + 2nc2H2 + 2nH2 = ^h- 


( 34 ) 
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These equations may be manipulated to obtain 
?t-ch 4 + ^co + 2nc2H2 
= nc (4ncH4 + 2nH20 + 2nc2H2 + 2) , 

«-H20 + nco = no (4ncH4 + 2nH20 + 2nc2H2 + 2). 

(35) 


(iBurrows & Shard Il999t iMoses et^ l2013ah . At 
low temperatures, methane do minates carbon monox¬ 
ide as the carrier of c arbon (iPrinn & Barshavl 119771 
[^^sha^&Xew^Q^^lthis tren^eyerses_ayii^i_tem- 
peratu res (iBurrows & Sharp! 119991 iLodders & FeglevI 

[ 200 l. 


Note that the number densities of the molecules marked by 
tildes have been normalised by nH 2 , while those of the atoms 
have been normalised by nn- The former are the mixing ra¬ 
tios, while the latter are the normalised elemental abundances. 

With two particle conservation equations and the expres¬ 
sions for K' and K 2 , we have four equations and four un¬ 
knowns. They can be manipulated to yield a cubic equation 
for the mixing ratio of methane, 

Co^ch4 + ('1'^ch 4 4" C2ficH4 + C 3 = 0; (36) 

which has the coefficients, 

Co = 2 K'K '2 (no - nc + 1), 

Cl —K' (4no — 4nc + 1) 

- K 2 [4nonc -f 2 (1 - 2no) (nc - 1)], (37) 

C 2 =2K' (no — n-c) — 4nc — 2no + 1, 

C 3 = — 2nc- 


While analytical solutions do exist for cubic equations, they 
possess multiple branches—some of which are complex— 
depending on tedious comb inations of the values of Co, Ci, 
C 2 and C 3 (iPress et al.ll2007h . Since it is difficult to determine 
a prior which solution branch ncH 4 is described by, we elect 
to solve the cubic equation using standard, canned numerical 
routines for solving polynomial equations. The other mixing 
ratios can be obtained via 


t^H 20 = 


2ho ( 7 C 2 ncH 4 + 2 ncH 4 + l) 


1 -I- Ar'ncH4 - 2no 
nco = 7 C'ncH 4%2 0 , 

^iC2H2 = 


(38) 


For completeness, we note that when acetylene is absent 
(762 = (()> the solution can be easily written down. 


«CH4 


-C 2 + (C| - dCiCa)'^' 
2Ci 


(39) 


Notice how the coefficient Co, and thus controls the ex¬ 
tent to which the mixing ratio of methane is described by a 
quadratic versus cubic equation. Physically, we expect that 
at low temperatures {K' <C 1), the mixing ratio of acety¬ 
lene is negligible. In this limit {K',K 2 ^ 1), we have 
ncHi ~ 2nc and ss 2no. These asymptotic solu¬ 

tions explain the relatively simple behavior of the mixing ra¬ 
tios at low temperatures, as seen in Figure |3] It also offers 
an easy explanation for methane and water switching roles as 
the dominant molecule when the carbon-to-oxygen ratio is ex - 
actly unity, as noted bv iKopparapu, Kasting & Zahnle|(l2M^ . 
since ncH 4 /«H 20 ~ nc/no. 

Figure|3shows the mixing ratios of methane, water, carbon 
monoxide and acetylene as functions of K'. Our analytical 
model produces the following, salient trends. 


• When the atmosphere has a solar abundance of ele¬ 
ments, water is always more abundant than methane 


• When the atmosphere is carbon-rich, methane is the 
dominant molecule but has to compete with acety¬ 
lene in some circumstances ( M adhusudh^ 120121: 
iMoses et akl l2013at IVenot et ^ 20151) . Water is the 
dominant oxygen carrier only at low temperatures, su- 
perceded by carbon monoxide at high temperatures 
(lMadhusudhanl2012tIMoses et al .1120 13 alibi) . 

• Cold atmospheres are always methane-rich at the ex- 
pense of carbon mo noxide, regardless of the C/O 
(iMadhusudhanI l2()T^ . The a bundance of water is es¬ 
sentially constant across C/O (lMadhusudhanl20 . 

• Hot atmospheres exhibit more complex behayior, in 
tha t they are methane-po o r and water-rich when C/O < 
1 (IMadhusudhanI 120121: IMoses et al.l IJoHH). For 
C/O > 1, they become methane-rich a nd water-poor 
(lMadhusudhani20l^IMoses et al.l2013a^ with methane 
dominating carbon monoxide as the carrier of car- 
bon when C/O becom es sufficiently larger than unity 
(IMadhusudhanI 120121) . When C/O is large enough, 
acetyle ne oyertakes methane as the dominant carrier of 
carbon (lMadhusudhanl20 1^ . 

These trends are in agreement with the nu merical calcu¬ 
lations of equilibrium chemistr y presented in IMadhusudhanI 
(l2(iTl) and IMoses et al.l (1201 3ah . but with one exception. We 
haye assumed K 2 IK' = 10 “^; higher yalues would produce 
the unphysical result that acetylene dominates carbon monox¬ 
ide, eyen at low temperatures (not shown). The cold model in 
the lower panel of Figure|3]shows an oyerabun dance of acety¬ 
lene, w hich is in disagreement with Figure 2 of iMadhusudhanI 
(l2012h . This discrepancy arises from the fact that K 2 IK' is 
not a constant and must possess a (steep) temperature depen¬ 
dence. 

To inyestigate this discrepancy further, we stop 
treating K' and as free parameters and instead 

relate them to temperature and pressure yia the 
Gibbs free energy tabulated in the JANAF database 
(http : //kinetics . nist. gov/ janaf /). Our imple¬ 
menta tion of this procedure is described in iHeng & LyonsI 
(I 2 OI 5 I) . In Figure 13 we recalculate the models in Figure 0 
We see that the basic trends preyiously discussed are pre- 
seryed, although the curyes display quantitatiye differences 
as expected. The mixing ratios yersus C/O match surprisingly 
well eyen at a quantitatiye ley el. The preyious result regarding 
acetylene is yerified to be an artifact of assuming K 2 IK' to 
be constant. We further yerified that acetylene becomes dom¬ 
inant oyer methane only for carbon-rich at mospheres with 
T > 2000 K (not shown), in agreement with IMadhusudhanI 

dmn . 

Oyerall, it is surprising how well our model is able to repro¬ 
duce the main trends of mixing ratios yersus the carbon-to- 
oxygen ratio. It is surprising because this rich yariety of be¬ 
hayior originates from the chemical analo gue of geometry— 
it is m erely stoichiometric book-keeping (ILodders & Fegle\^ 
l2002h . The dependence of the normalised equilibrium con- 
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slants on temperature is a distraction if all one seeks is to un¬ 
derstand these trends in a qualitative sense. 

5. SUMMARY 

We have presented a unified, novel, self-consistent formal¬ 
ism for understanding the atmospheric chemistry of exoplan¬ 
ets from the viewpoint of an astrophysicist. In doing so, we 
addressed ambiguities associated with the equilibrium con¬ 


stant and obtained a novel derivation of the Arrhenius equa¬ 
tion. We also generalised previous work on analytical mod¬ 
els of systems in chemical equilibrium with carbon, hydrogen 
and oxygen and showed that they reproduce several key trends 
published in the literature and computed using more sophisti¬ 
cated numerical calculations. We anticipate that such models 
are useful for inclusion in retrieval models of exoplanetary 
atmosphe res to maintain their chemical plausibility as a first 
approach (lBennekell2015h . 
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APPENDIX 

CHEMICAL KINETICS 

For completeness, we restate the formalism concerning chemical kinetics. 


Evolution equations 

Unlike Gibbs free energy minimisation, chemical kinetics is the treatment of a network of reactions as a system of mass 
conservation equations. The evolution of the reaction Xi is described by the partial differential equation. 


1 (dnxi ^ 

ai V dt "" dx^ ) 




Jxi- 


(Al) 


The importance of the 1/ai factor cannot be over-stated: it allows for the reaction rates of reactants and produ cts with different 
stoichiometric coefficients to be placed on the same footing (lJohnstonlll966t [Steinfeld, Francisco & Haseiri989h . The production 
and loss rates are 


V = 


(A2) 


The reaction rate associated with photochemistry is given by jTxi. which generally depends on nxi ■ 

The diffusion coefficient (K^z) is used to mimic advection, convection and turbulence and subsume their collective influence 
into a single free parameter. Generally, advection, convection and turbulence hardly resemble diffusion in any rigorous way—one 
often argues that these processes operate on scales that are so small, compared to the characteristic atmospheric length scale of 
interest, that it “looks” like diffusion. The use of Kzz is rigorous and exact only for molecular diffusion. Notwithstanding, the 
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inclusion of a diffusion coefficient allows us to treat situations with disequilibrium chemistry induced by atmospheric motion or 
mixing without resorting to a full-blown, three-dimensional calculation. 

For the product Z\, the evolution equation is 


1 (d'^nz^ \ 
bi V dt "" dx^ ) 


= V' - C'nl 




where the production and loss rates are 


V — 
C = 


(A3) 


(A4) 


Why photochemistry is a disequilibrium effect 

In the absence of atmospheric mixing (Kzz = 0), we may add the evolution equations for Xi and Zi to obtain 

1 9nxi , 1 dnzi 


-f 


= —fTxi — Jzi ■ 


tti dt bi dt 

If we integrate this expression, we obtain 

nxi , nzi 

Oi 

If we do the same for all combinations of reactants and products, then we obtain 


H —— — J (i7xi + <Jzi) dt + C. 


nxi ^ 

tti 02 bi 62 


— ~ J {Jxi + J7x2 + dt -\- C , 


(A5) 


(A6) 


(A7) 


where C and C are constants of integration. 

This result informs us that photochemistry is an intrinsically disequilibrium effect, because it allows the total number of 
particles in the system to vary with time. In its absence, the total number of particles is an invariant quantity. 


Producing chemical equilibrium in the steady-state limit 
If we neglect atmospheric mixing and photochemistry, the steady-state limit of the evolution equations yields 




(A8) 


Since this is identical to the setup in which we used to define our dimensional equilibrium constant we conclude that our 

evolution equations correctly produce chemical equilibrium in the steady-state limit. 


GIBBS FREE ENERGY FOR HYDROGEN ATOM 


We use the Gibbs free energy associated with the hydrogen atom from the JANAF database. Here, we list it in units of kJ/mol/K, 
from 0 to 6000 K (in increments of 100 K) and at Po = 1 bar; 216.035, 212.450, 208.004, 203.186, 198.150, 192.957, 187.640, 
182.220, 176.713, 171.132, 165.485, 159.782, 154.028, 148.230, 142.394, 136.522, 130.620, 124.689, 118.734, 112.757, 
106.760, 100.744, 94.712 , 88.664, 82.603, 76.530, 70.444, 64.349, 58.243, 52.129, 46.007, 39.877, 33.741, 27.598, 21.449, 
15.295, 9.136, 2.973, -3.195, -9.366, -15.541, -21.718, -27.899, -34.082, -40.267, -46.454, -52.643, -58.834, -65.025, -71.218, 
-77.412, -83.606, -89.801, -95.997, -102.192, -108.389, -114.584, -120.780, -126.976, -133.172, -139.368. If we denote each of 
these numbers by Gh, then we have AGq = —2G'h for the net reaction in equation (l24l i. It follows that 


K = — exp —; 


Pn 


7? • T 


(Bl) 


where T^univ = 8.3144621 J K ^ mol ^ is the universal gas constant. See iHeng & Lvon^ (120151) for more explanation on the 
unit conversion between AGq and AGq- 











